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Statistical properties of non-symmetric real random matrices of size M, obtained as truncations of random 
orthogonal TV x TV matrices are investigated. We derive an exact formula for the density of eigenvalues which 
consists of two components: finite fraction of eigenvalues are real, while the remaining part of the spectrum is 
located inside the unit disk symmetrically with respect to the real axis. In the case of strong non-orthogonality, 
M/N =const, the behavior typical to real Ginibre ensemble is found. In the case M = TV — L with fixed L, a 
universal distribution of resonance widths is recovered. 
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Random unitary matrices [1] are used in numerous phys- 
ical applications including chaotic scattering, conductance in 
mesoscopic systems Pf] or statistical properties of periodically 
driven quantum systems y[]. In several applications one needs 
to restrict this class to real, orthogonal matrices, which are as- 
sumed to be distributed uniformly with respect to the Haar 
measure on the orthogonal group. This is the case while de- 
scribing quasiparticle excitations in metals and superconduc- 
tors J4j,|5|], or quantum maps performed on real quantum states 
|g]. In applications one is often led to consider square trun- 
cations of random unitary matrices of large matrix dimension 
TV. These matrices have been used to describe quantum sys- 
tems with absorbing boundaries [7] and found applications in 
various physical problems including optical and semiconduc- 
tor superlattices |8], problems of quantum conductance |91, 
distribution of resonances for open quantum maps H 1 OM 1 311 . 
The operation of truncation does not preserve unitarity and as 
a consequence the eigenvalues move inside the unit disk in the 
complex plane, see, e.g., HQ, and the emerging eigenvalue 
statistics appear to be universal Holl . 

The aim of this work is to present a comprehensive study 
of truncations of random orthogonal matrices. We consider 
an ensemble of TV x TV real orthogonal matrices Q with flat 
matrix distribution on the orthogonal group. The goal of this 
paper is to compute analytically the full distribution of eigen- 
values of the top M x M sub-matrix A of Q = j|) for 
any values of M and TV II15ll . This then gives the eigenvalue 
correlation functions in a closed form and we investigate two 
asymptotic regimes of direct physical importance. For sim- 
plicity we consider M even, however, our results can be gen- 
eralized to odd dimensions as well. 

The orthogonality condition for the first M columns of Q, 
A T A + G T C — 1, defines a manifold in the space of real 
TV x M matrices which can be identified with the flag mani- 
fold 0(N)/0(L) where L = TV - M. Correspondingly, the 
density of distribution of A, C is 



P(A,C) 



vol O(L) 



8{A T A + C T C - 1) 



volO(TV) 

where volO(TV) = 2 N njli ^ /2 /r(j/2) is the volume of 



the orthogonal group O(N) |@]. The density function of A is 
P(A, C) integrated over C. For L > M the C-integral can be 
performed by scaling C — > Cy/l — A T A. The corresponding 
Jacobian is | dct(l — A T A)\ L / 2 and from the <5-function one 
gets another factor | dct(l — A T A)\~( M+1 ^ 2 . The remaining 
integral yields the volume of O(L) fO(L — M), leading to[9] 



P(A) 



vol 0(L) vol O(L) 
volO(TV) vol O(L-M) 



det(l-A T AY 



(1) 



If L < M the probability density function of A contains sin- 
gular (/5-function) terms due to finite mass of the boundary of 
the matrix ball A T A < 1 which supports the distribution of A. 

By setting M = 1 in Eq. (Q~|) one recovers the distribution 
P(A) oc {I- A 2 )— of a single entry of Q. Obviously, for 
TV large P(A) can be approximated by a Gaussian. This is 
also true for matrix blocks if the size of truncation M <C TV. 
Eigenvalue correlations in such a Gaussian regime are now 
accessible following the recent progress [[16^2011 made for the 
real Ginibre ensemble ll2lll . In this paper we investigate two 
non-Gaussian regimes: (i) strong non-orthogonality M, L oc 
TV; and (ii) weak non-orthogonality M oc TV and L <C TV. In 
the former we again recover Ginibre eigenvalue correlations 
while the latter yields eigenvalue statistics from a different 
universality class 

We shall obtain the joint distribution of eigenvalues of A 
directly from P(A, C) reducing the general case to M = 2. 
Thus, it is instructive to consider first Al = 2. In this case, one 
can bring A by a rotation to the form A = {_g\^) with real 
Ai,2 and 5 from which the eigenvalues Zi of A can be recov- 
ered via the relations det A — Ai A2 + 8 2 and Tr A = Ai + A2 
1 1811 . The matrix A can have either two real eigenvalues 
in which case we choose zi > z 2 or two complex conju- 
gate eigenvalues in which case we choose Imzi > Imz 2 . 
The joint distribution d/^(z 1; z 2 ) for L > 2 in each of these 
two cases can be obtained from Eq. (Q]) by making use of 
det(l - A T A) = (1 - zf)(l - z 2 ) - AS 2 and integrating 
out 8 subject to the constraints 4<5 2 < (1 — zf)(l — z 2 ) and 
AS 2 > — (z\ — Z2) 2 arising due to the chosen parameterization 
of A, see Hill , and the positivity of 1 — A T A. The resulting 



2 



expression can be presented in a form that encodes both cases 
(real and complex): 

dfj,(zi,z 2 ) = (zi - z 2 ) f{zi)f{z 2 )dz\ A dz 2 (2) 



with (z = x + iy) 

f{z) = HL_h l{ ___ . 2iL . 2 

Zir 



du(i- u 2 y 



(3) 



If the eigenvalues of A are real then the integral in (0 yields 
a Beta function. By integrating dfi(zi,z 2 ) over the triangle 
— \<z 2 <z\< \ one gets the probability for A to have 
two real eigenvalues: p% L = B(L/2 + l/2,£)/B(i/2,£ + 
1/2). In the limit i — > oo this converges to l/\/2 which, as 
expected, is the value for the real Ginibre ensemble. 

One recovers d/x(zi, z 2 ) for L = 1 by letting L — » 1 in (0. 



In this limit / 2 (z) = (2tt|1 



As a simple check one 



can evaluate the probability v 2 l=\ f° r A t0 have two complex 
eigenvalues. In this case z\ = z 2 * so that (z% — z 2 )dz\l\dz 2 = 
Aydxdy in (|2). Integrating d//(zi, z 2 ) over the upper half of 
the unit disk \z\ < 1 one gets p§ 1 = 1 — 2/7r = 1 — p^, in 
agreement with the above result for v 2 l- 

We now proceed with the general case of M even. By an or- 
thogonal transformation O drawn fromO(M)/0(2) M / 2 one 
can bring A to upper triangular form A — 0(A + A)O t , 
where A is block diagonal with 2x2 blocks Aj in the diago- 
nal and A is block triangular with 2x2 blocks Ajj (I < J) 
above the diagonal. Integrating P(A, C) over O, A and C we 
obtain the distribution of A, 

P(A) = (47r)- M / 2 CMxA'(A)x (4) 

J dAj dCS((A T + A T )(A + A) + C T C - 1) 

mthc M ,L=volO(L)volO[M)/volO(L+M). Here A'(A) = 
A(A)/J}j A(Aj) is the Jacobian of the transformation from 
A to O, A and A, see J3, and A(A) = U l<:j (zi - Zj) is the 
Vandermonde determinant of the eigenvalues Zj of A. Thus 
A' does not count the pairs i < j from the diagonal blocks. 

Writing C as C = [Ci, . . . , Cm/ 2 ] in terms of L x 2 sub- 
blocks Cj, we have for I < J the equations 



AJA 
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K<I 



AkiA kj — —CfCj 



(5) 



due to the ^-function in (01), which start with AfAu = 
-CjCj for J > 2, A|A 2J + Aj^j = -CgCj for J > 3, 
etc. Because of the triangular structure of A these equations 
can iteratively be solved for Ajj, Thus integration over Au 
can be performed with a Jacobian Yli< j I det A-H 2 - There 
remains integration over Cj\ 



M/2 

l\[dCj 5(A T jAj + 



J=l 



K<I 



KI A KJ + CjCj-l)]. 



From Eqs. (|5]l one sees that A/j for I < J is linear in Cj 
and depends otherwise on Ck for K < J. Thus we have an 



integral / iffli dCj 6{A T jAj + CjXjCj - 1) where Xj is 
a linear operator acting on Cj per matrix multiplication. The 
scaling Cj — > (1/s/Xj)Cj results in a Jacobian that cancels 
the Jacobian Il/< j I det A/|~ 2 from the previous integration 
over Au and we end up with 



P(A) _ A'(A) 



cm.l (4tt)^/ 2 



.M/2 

/ HdCjS(A^Aj + CjCj-l). 



where Aj's are 2 x 2 and Cj's are L x 2. Now we proceed 
calculating P(A) using our result for M = 2. Performing 
integration over Cj yields 



P(A) A'(A) 



call (4tt) m /2 11 volO(L-2) 



M/2 

n 



vol O(L) 



det(l-AjAj)^, 



and we obtain the joint (full) distribution of eigenvalues of A: 

M 

dfi(z 1 ,...,z M ) = Z M ] L Yl {Zi - Zj) Y[f(Zj)dzj (6) 



l<i<j<M 



with Z M | L = [ J L!/(27r) i ] M / 2 CAf,L, dzj = cfoj + idy^ and 
ordering Z\ > z 2 > ... > Zm if all eigenvalues are real, 
zi = z 2 *, Imzi > 0, Z3 > Z4 > . . . > zm if two eigenvalues 
are complex conjugate and so on. Eqs. (0, (O represent our 
first main result. 

Eq. © is of the same form as the corresponding one in the 
real Ginibre ensemble II 1 811 but with a different weight func- 
tion /(z). Following [18] one can obtain the eigenvalue cor- 
relation densities in terms of a skew-symmetric kernel 



M M 



K(z 1 ,z 2 ) = J222( A ~ 1 ^ z r 



k=l 1=1 



fe - 1 - i r 1 withd 2 ; 



where Am = J Jd 2 z 1 d 2 z 2 J : (z 1 ,z 2 )z 1 
dxdy and 

F{zi,z 2 ) = /(zi)/(z 2 )x 

[2i6 (2) (zi - z 2 *)sgn(yi) + 5{y l )5{y 2 )sgn{x 2 - a?i)], 

with S^(zi — z 2 *) — 8{x\ — x 2 )S(yi + y 2 ). For example, 
the one-point density is given by an integral 



Ri{z) 



d 2 z 2 F(z, z 2 )K[z 2 , z). 



(7) 



Higher order densities R n are given by a Pfaffian involving 
T and JC 111 811 . Expressing the kernel JC in terms of skew- 
orthogonal polynomials with respect to weight function T 
leads to an alternative Pfaffian representation for R n B16lll7ll . 

Evaluating the kernel K, (or, equivalently l24ll . the corre- 
sponding skew-orthogonal polynomials) in a closed form is 
an important step on the way to obtaining eigenvalue statistics 
ifnll . In our case the kernel can be found by exploiting its re- 



lation to averages of the characteristic polynomials IU8LI25H of 
truncations of orthogonal matrices of smaller dimension: 

fc(zi,z 2 ) = (zi - z 2 ) (det(zi - A)(z 2 - A)) A 
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where the matrices A are square truncations of size M — 2 of 
random orthogonal matrices of size TV — 2. Due to the invari- 
ance of the distribution of A the above average is effectively 
an average over the eigenvalues of A T A: 



(det( Zl -A)(z2-A)) A =^2( 



2122) 



m=0 



\ c n— m 

(A T A)) A 

^n—m (1) 



where n = M — 2 is the size of A and e m (X) is the m-th 
elementary symmetric function in eigenvalues of X. The in- 
tegral [e m {A T A)) A can be reduced 12611 to Selberg's integral 
02911 . yielding the kernel in a closed form, 



£(2l,2 2 ) = ( 2 1 - z 2) 



if- 

m=0 



m)! 



L! m! 



(2l2 2 ) 



(8) 



which is our second main result. The truncated binomial se- 
ries on the rhs can be expressed [23] in terms of an incomplete 
beta function which comes in handy for asymptotic analysis of 
eigenvalue statistics in the limit of large M. 

The one-point eigenvalue density (0 is composed of two 
parts: Ri(z) = p c {z) + S(y)p R (x), where p c (z) = 
2f{z) 2 \JC(z, z*)\ is the density of complex eigenvalues and 



dx 2 sgn(x 2 - x)IC(x 2 ,x)f(x 2 )f(x) (9) 



is the density of real eigenvalues of truncated orthogonal ma- 
trices, with the normalization f\ x \ <x Ri{z) d 2 z = M. 

We shall first look at real eigenvalues of A. Integrating 
Pr(x) over — 1 < x < 1 one obtains the expected total num- 
ber Nr of real eigenvalues, 



N R = 1 



ds 



T I S 2 (L/2,1/2) (i + l,M-l), 



where I x (a, b) = j* - t) h ~ l dt/B(a, b) is an incom- 

plete Beta function. For large matrix dimensions (M ^> 1) 
Nr, in the leading order, is described by a simpler expression: 



N R 



2 artanh y/M/N 
B(L/2,l/2) 



Depending on L this leads to different scaling laws for Nr. In 
the limit of strong non-orthogonality when both M, L oc N, 
the number of real eigenvalues oc y/M which is characteristic 
of the real Ginibre ensemble 127]. On the other hand, in the 
limit of weak non-orthogonality when L -C M, Nr grows 
logarithmically with M, N R ~ (logM)/B(L/2, 1/2). 
Eq. (O can be transformed to a form, 



Pr{x) 



S(£/2,l/2) l-x 2 
B(M/2,L/2) 



7 x2 ((M-l)/2,(L + 2)/2), 




FIG. 1. Rescaled density of real eigenvalues p — for M = 32 
and L = 2, 4, 8, 16. The inset shows the number of real eigenvalues 
Nn for M — 32 as a function of L obtained from samples of 10 3 
realizations (squares) and analytical prediction (solid line) 



revealing the Artanh Law of the distribution of real eigenval- 
ues for large matrix dimensions M (p = M/N): 



1 



1 



B(L/2,l/2) 1 



y/Jl<X < y/JI. 



In the limit of strong non-orthogonality p < 1 and pr(x) 
is supported strictly inside the interval (—1, 1), vanishing at 
a Gaussian rate at the boundaries x = ±«//I. In the limit of 
weak non-orthogonality p 1 and pr(x) becomes singular 
at the boundaries x — ±1 of the real eigenvalue support. Cor- 
respondingly, real eigenvalues accumulate close to x = ±1. 
Setting x = 1 — u/M and taking the limit M —> oo one ob- 
tains the density profile near the accumulation point x = 1: 
Pr{x)/M ~ p(u) with 



p(u) 



2T(L/2) 
(2u) L 



(1-ut+V(£/2 + 1,u)) + 



- 7 *(L + 1,2 U ), 



B(L/2,l/2) 

where 7* (n, x) = (x™r(n)) -1 t™ _1 e~*di is an incomplete 
Gamma function. For small u, p(u) ~ (2r(-|)) _1 ii"5" _1 , 
showing a transition, as L increases, in nature of the bound- 
ary points x — ±1 from 'attractive' (L — 1) to 'repul- 
sive' (L > 4), see Fig. [T] On the other hand, for large u, 
p(u) ~ (B(L/2, l/2)u)~ 1 , exhibiting a heavy tail which 
manifest itself in the log M asymptotics for Nr. 

Now we shall look at complex eigenvalues. Their density 
Pc(z), visualized in Fig.|2j vanishes on the real line. In the 
limit of strong non-orthogonality, close to the real line the 
complex density is described by the same scaling law 



p c (z)~pR(x) 2 h(ypR(x)), h(y) = 4n\y\e^y erfc(V4^|y|) 



as in the real Ginibre ensemble 1 28], while away from the real 
axis pc (z) is the same as for truncations of random unitary 



matrices 112211 . pc(z) 



z\ 2 ). In the limit 



of weak non-orthogonality away from the real axis 

pc(z) ~ v 2 h(A-nv{\ - \z\)), h(u) = ATtLu L - 1 1 *{L + l,u), 
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FIG. 2. Contour plot of the normalized density of complex eigen- 
values pc(z)/(l — Nr) of truncations of random orthogonal ma- 
trices for a) TV = 12 and M = 10 (weak non-orthogonality) and 
b) TV = 12 and M — 6 (strong nonorthogonality) with the larger 
values coded darker. The density plotted in the upper half of the 
complex planes represents analytical results, while the density in the 
lower part is obtained numerically for a sample of 0.5 x 10 r random 
orthogonal matrices. 



where v 



M ■ 



is the density of distribution of eigenvalues of 



random orthogonal matrices along the unit circle. Again, the 
limiting density profile does not depend on the angle and is 
the same as for truncations of random unitary matrices [22]. 

Eigenvalue correlations for truncations of random or- 
thogonal matrices can also be obtained in a closed form. 
In the limit of strong non-orthogonality, after appropri- 
ate rescaling the eigenvalue correlations become identical 
to those in the real Ginibre ensemble. For example, at 
the origin, f 2 (z/y/I) ~ ^/ L/^-k)^-^ erfc(V%|) and 
K(zi/VL, Z2/VL) ~ ^7j(zi — Z2)e zlZ2 which is what one 
gets in the real Ginibre ensemble lfl7r - [l9ll . The limit of weak 
non-orthogonality away from the real axis the eigenvalue cor- 
relations for truncations of random orthogonal matrices are 
exactly the same as those found for truncated random unitary 
matrices B22l 12311 . New correlation laws arise in the vicinity 
of x = ±1. 

TABLE I. Ensembles of random nonhermitian matrices 

matrices complex real /j — M/N 



Haar 
measure 

truncations of 
matrices 

Ginibre 
ensemble 



U(N) O(N) fi = 1 
unitary orthogonal 1 > fx > 
complex real /j, — } 



In conclusion, we have found the full probability distribu- 
tion of eigenvalues of truncated random orthogonal matrices 
and obtained the eigenvalue density and higher order corre- 
lation functions in a closed form. This work completes our 
understanding of ensembles of non-hermitian random matri- 
ces which are derived from random unitary and orthogonal 
matrices. As shown in Table I truncations of random orthog- 



onal (unitary) matrices form an ensemble which interpolates 
between matrices distributed according to the Haar measure 
on the orthogonal (unitary) group 1E9I1 and the real (complex) 
Ginibre ensemble 112 111 . 
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Note added. After submitting our work we learned about a 
very recent preprint of Forrester ll3Qfl extending our study 
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